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Abstract: Experiments at electron-positron colliders can search for dark matter par- 
ticle pair-production in association with a photon. We estimate the sensitivity of 
this search at the proposed International Linear Collider (ILC), under a variety of 
run scenarios. We employ the effective operator formalism to provide a quasi-model- 
independent theoretical description of the signal, and present the reach of the ILC in 
terms of the scale A suppressing the dark matter-electron coupling operator. We find 
that at the 250 GeV center-of-mass energy, the ILC can probe A up to 1 — 1.2 TeV, a 
factor of 2.5 — 3 above the best current bounds from LEP-2. With 1 TeV energy and 
polarized beams, the reach can be extended to 3 — 4 TeV. The ILC can discover this 
signature even if annihilation to electrons provides only a small fraction of the total 
dark matter annihilation rate in the early universe. We also argue that large regions of 
parameter space allowed by current LHC and direct detection bounds will be accessible 
at the ILC. 
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1 Introduction 



While the existence of dark matter has been firmly established through observation of 
its gravitational effects on multiple scales, its microscopic nature remains unknown. It 
is clear that none of the Standard Model (SM) particles can account for dark matter, 
so that new physics must be involved. This situation motivated a variety of searches 
for dark matter through non-gravitational interactions, such as direct, indirect and 
collider signatures. No firm detection has been achieved as yet, but several intriguing 
hints have been reported. 

The focus of this paper is on collider signatures of dark matter. From the theo- 
retical point of view, a very attractive class of dark matter candidates is the weakly 
interacting massive particles, or WIMPs. WIMPs fit naturally in a variety of popu- 
lar extensions of the SM at the electroweak scale, and they generally have the right 
order-of-magnitude thermal relic density to fit the observed dark matter abundance 
(the so-called "WIMP miracle"). It is plausible that WIMPs could be produced at 
significant rates at current or near-future high-energy colliders, giving potentially ob- 
servable experimental signatures. Extensive studies of such signatures appeared in the 
literature. However, most studies are performed in the context of a particular particle 
physics model, such as supersymmetry (usually with very specific superpartner spec- 
trum in mind), Little Higgs, etc. Given that a large number of particle physics models 
contain WIMPs, it is highly desirable to identify and focus on signatures which do not 
depend sensitively on model details, but are as model-independent as possible. 

The first example of a model-independent approach to collider signatures of dark 
matter appeared in 2004 in Ref. [1]. This paper argued that the process of pair-creation 
of dark matter particles in, for example, e + e~ collisions, is just the reverse of one of 
the annihilation processes which determine the thermal relic density of the WIMPs. 
This allows one to use the measured abundance of dark matter to predict the rate of 
its production at an e + e~ collider. An observable signature is obtained if a photon 
is emitted before the collision: since the WIMPs escape the detector unobserved, the 
signature is 7 + f! . Using soft /collinear factorization theorems, the rate of this process 
can also be predicted in a model-independent way, up to a small number of unknown 
phenomenological parameters, each with a clear physical significance. Detailed Monte 
Carlo studies of the sensitivity of the proposed International Linear Collider (LLC) to 
this process have since been performed [2, 3]. 

More recently, an alternative theoretical framework for a quasi-model-independent 
studies of dark matter at colliders was proposed [4]. The starting point is a La- 
grangian in which the dark matter is coupled to the SM through a set of effective 
(non-renormalizable) operators. Such operators represent the effect of heavy "media- 
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tor" particles, but all dependence on the specifics of the mediators is contained in a 
small number of operator coefficients. Most of the applications of this approach have 
focused on dark matter signatures at hadron colliders, j + fir and 7 + ^. These are in- 
duced by the same operators as the dark matter signal in direct detection experiments, 
so correlations between the two types of experiments can be understood. Tevatron 
and LHC experiments have already reported bounds from searches performed with this 
approach [5] . Effective operator formalism can also be applied to dark matter produc- 
tion in e + e~ collisions. This has been done in Ref. [6], which reported bounds on dark 
matter-SM couplings from non-observation of anomalous 7 + f events at LEP-2. 

An energy frontier e + e~ collider, the International Linear Collider (ILC), is one 
of the most compelling proposals for future major facilities in fundamental physics. 
The most immediate physics motivation for such a machine is to study the recently 
discovered 125 GeV boson in detail, and confirm or disprove its identification with the 
SM Higgs. For this task, the center-of-mass energy of e + e~ collisions of about 250 GeV 
would be ideal. In the future, the accelerator can be upgraded to reach the collision 
energies of 1 — 1.5 TeV. One of the most intriguing experiments that can be performed 
at the ILC is a search for dark matter. Our goal in this paper is to evaluate the potential 
sensitivity of such a search, under a variety of machine parameters suggested by recent 
design studies. To minimize dependence on unknown particle physics, we will employ 
the effective operator formalism mentioned above. 

Given the time scale for the completion of design and construction of the next e + e~ 
collider, it is quite possible that WIMP dark matter will be discovered, for example 
at the LHC or via direct detection, before the experiments discussed in this paper can 
begin. In this scenario, the main task of the e + e~ collider will be to perform detailed 
measurements of the dark matter particle properties. Still, the dark matter signal will 
need to be established first, so our study would still be relevant. Moreover, it has 
been seen in similar examples that the reach for discovery and model discrimination 
is typically roughly similar (see e.g. [7, 8]), so our study gives at least a rough sense 
of how well different models (for example, different effective operators for SM-WIMP 
coupling) can be distinguished. 

The rest of the paper is organized as follows: We describe the theoretical framework 
of the analysis, the collider run scenarios considered, and other assumptions, in Sec. 2. 
Our estimates of the dark matter search reach at the ILC are presented and discussed 
in Sec. 3. We summarize the main conclusions of the analysis in Sec. 4. Appendix A 
contains the analytic formulas for the differential cross sections of the radiative WIMP 
production in (polarized) e + e~ collisions, which have not previously appeared in the 
literature. 
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Figure 1. Feynman diagrams for radiative WIMP pair-production in e + e collisions, in the 
operator formalism. 



2 Setup 

While the operator formalism can be used for WIMPs of any spin, we will assume, 
for concreteness, that the WIMP is a spin-1/2, Dirac fermion \. The coupling of the 
WIMPs to electrons and positrons has the form 

£tat = ^0i, (2.1) 

where A roughly corresponds to the energy scale of new physics that provides the 
coupling, and Oi is one of the following four- fermion operators [6]: 

O v = {xiMh^), (vector) 

Os = (xx) (P?) ) (scalar, s — channel) 

O a = (xi^l5X)(h^ 5 £) , (axial - vector) 

O t = (xt)$x), (scalar, t-channel). (2.2) 

The notation in parenthesis describes the simplest kind of a mediator particle that 
would induce each operator. We will always consider the case when the mediator mass 
is well above the collision energy y/s, and our results will not depend on how the opera- 
tors (2.2) are induced; the names are only used as a convenient way to label operators. 
Since the WIMPs do not interact in the detector, the 2 — > 2 process e + e~ — > XX is 
invisible; an extra "tag" particle needs to be added to the final state to make it observ- 
able. A photon can always be emitted from the initial state independently of the nature 
of the WIMPs and their couplings, making it a robust choice for the tag particle [1]. 
We will thus consider the process e + e~ — > XXli mediated by Feynman diagrams in 
Fig. 1, and leading to the observable 7 + f, final state. We have computed the double- 
differential cross sections, dE d ^ osd , analytically for each of the four interactions listed 
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ILC-250 


ILC-500 


ILC-1000 


ILC-500P 


ILC-1000P 


y/s, GeV 


250 


500 


1000 


500 


1000 




250 


500 


1000 


250 


500 


P-/P+ 


0/0 


0/0 


0/0 


+0.8/+0.5 


+0.8/+0.5 


cr b g, fb 


1.07 


1.74 


2.36 


0.21 


0.25 



Table 1. Accelerator parameters for each of the run scenarios studied in this paper, and the 
SM background cross sections within the experimentally accessible region. 



in (2.2) and for all possible combinations of electron and positron beam polarizations. 
The formulas are presented in Appendix A. 

The main irreducible background to the search for the 7 + f) signature is the 
SM process e + e~ — > vv^f. We have evaluated the cross section for this process using 
MadGraph/MadEvent v5 .0 [9] . We ignore instrumental backgrounds, and assume a 
systematic error on the background prediction of 0.3%. We combine the statistical and 
systematic errors in quadrutures to obtain the significance: 

Sig = Nsi& =. (2.3) 

y/N hg + (0.003iV bg ) 2 

The reach of a given experiment is estimated by requiring Sig> 3. In all scenarios 
considered below, Ab g 3> 1, so the use of Gaussian statistics is well justified. We present 
results for several sets of collider parameters, summarized in Table 1. For the proposed 
ILC, we follow [10] and assume three stages of operation: ILC-250 (y/s = 250 GeV, 
L int = 250 fb" 1 ); ILC-500 {y/s = 500 GeV, L- mt = 500 fb" 1 ); and ILC-1000 (y/s = 1 
TeV, L int = 1000 fb -1 ). We assume unpolarized beams in these three scenarios. In 
addition, we consider the possibility of beam polarization at 500 GeV and 1 TeV, taking 
the electron and positron polarizations to be P_ = +0.8 and P + = +0.5, respectively. 1 
(Note that in the polarized beam case, the uncertainty in the degree of polarization 
should be included in the systematic error; the 0.3% error assumed in our study would 
require knowledge of beam polarization at this level or better.) On the detector side, 
the parameters relevant for our study are the minimum photon energy L™ m , and the 
minimum value of the angle of the photon to the beamline, characterized by | cos # 7 | max , 
for which the photon can be detected. For all scenarios, we use E™ m = 8 GeV and 
I cos^l 111 ^ = 0.995. In addition, to further suppress the background, we discard events 
with photon energies corresponding to recoil against an on-shell Z. Specifically, we 



^^Our notation for polarization is discussed in Appendix A. The signs are chosen so that electrons are 
predominantly right-handed and positrons left-handed, suppressing the contribution of the diagram 
with the t-channel W exchange to the background. 
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define z = 2£ 7 /^, and discard events with z E [0.8,0.9] for ILC-250, z e [0.95,0.98] 
for ILC-500(P), and z G [0.98,0.99] for ILC-IOOO(P). The cross sections of the SM 
background process e + e~ — > Z/P7, after all cuts, are listed in Table 1. 

3 Results 

Our estimates of the reach of the ILC dark matter searches, with parameters described 
above, are given in Fig. 2. In terms of the scale A, the accessible regions are approxi- 
mately independent of the dark matter particle mass, until the kinematic reach of the 
collider (roughly M x » yfs/2) is reached. The dependence on the helicity structure of 
the operator is significant but not overwhelming, with all bounds within roughly 30% 
of each other. A 250 GeV ILC (the "Higgs factory") is sensitive to scales A up to about 
1 — 1.2 TeV, a factor of 2.5 — 3 higher than the LEP bounds given in Ref. [6]. Fur- 
ther improvements are possible by increasing collision energy and/or using polarized 
beams. For example, a 1 TeV collider with polarized beams can probe A up to 3 — 4 
TeV, depending on the operator. In all cases, the reach is significantly higher than 
the center-of-mass energy in e + e~ collisions, justifying the use of the effective operator 
formalism in our analysis. 

A lower bound on A is equivalent to an upper bound on the annihilation cross 
section XX e + e~. Let us assume that x is a thermal relic, and that e + e~ is the 
dominant annihilation channel in the early universe. The thermally-averaged value 
of this cross section, evaluated at the temperature at which decoupling takes place, 
determines the present density of WIMPs. In Fig. 4, we plot the projected ILC upper 
bounds on the quantity (cr(xx ~~ e + e~)t>), where v is the relative velocity of the 
colliding WIMPs and the brackets denote thermal averaging. (We assume (v 2 ) = 0.24 
at freeze-out.) The observed present density of WIMPs requires (<r(xx ~~ >* e + e^)v) m 
3x 10 -26 cm 3 /sec. It is clear that the ILC will be able to probe significantly smaller cross 
sections, throughout the kinematically accessible WIMP mass range. In particular, in 
the low WIMP mass region, between 1 and 10 GeV, the ILC-250 will be able to discover 
the WIMPs even if their annihilation to electrons provides only ~ 10~ 5 — 10~ 3 of the 
critical value. The LEP bounds on the annihilation cross section [6] would be improved 
by one to two orders of magnitude. 

It is also interesting to compare the ILC reach with that of direct detection searches 
for dark matter, as well as searches at the LHC. The LHC searches for WIMP pair- 
production (in association with a jet or a photon) in pp collisions. Direct detection 
experiments are primarily sensitive to WIMP scattering on nuclei. Both signals are 
dominated by WIMP-quark interactions, so any comparison with the ILC bounds re- 
quires making theoretical assumptions regarding the relation between WIMP-quark and 
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Figure 2. Reach of the ILC dark matter searches, for the 5 run scenarios defined in Table 1, 
as a function of the dark matter mass M x and the scale A of the WIMP-electron coupling 
operators. The regions below the curves are accessible at the 3-sigma level. 



WIMP-electron couplings. Two distinct scenarios are plausible. First, the WIMPs may 




Figure 3. The one-loop diagram inducing the WIMP-quark coupling from a tree- level WIMP- 
electron operator. 



couple to quarks directly, at tree level, through operators similar to those in Eq. (2.2). 
Second, the WIMP may only couple directly to leptons, in which case the dominant 
coupling to quarks appears via the one-loop diagram shown in Fig. 3. (Which of these 
scenarios is realized depends on the quantum numbers and couplings of the new parti- 
cles at scale A that induce the WIMP couplings to the SM.) Furthermore, while collider 
bounds are relatively insensitive to the structure of the operator, the bounds from di- 
rect detection vary widely depending on this structure. Let us discuss a few interesting 
possibilities. 

Direct WIMP-Quark Couplings: If the WIMP-quark operator mediates spin- 
independent WIMP-nucleon scattering at low velocities, such as scalar or vector oper- 
ators, then the bounds from direct detection are very strong. The current XENON100 
bound [11] rules out operators of this kind up to A ~ 20 TeV. This is well above the 
ILC sensitivity even in the most optimistic scenario, but it is worth emphasizing two 
points. First, XENON100 (and most other direct detection experiments) rapidly lose 
sensitivity for WIMP masses below 10 GeV, and, as is already well known, collider 
searches provide the more sensitive probe in that region. Second, there is no a priori 
reason why the scale suppressing the quark and lepton operators should be the same. 

If the WIMP-quark operator only mediates spin- dependent WIMP-nucleon scat- 
tering (e.g. axial- vector or tensor couplings), then the direct detection bounds on A 
are much weaker, of order a few hundred GeV. ILC searches will be sensitive to much 
higher scales, and in addition will cover the region of low-mass WIMPs. 

For all operator structures, the current LHC bounds on A are of order 700 GeV, 
well below the ILC sensitivities. (Although, once again, we emphasize that there is no 
a priori reason why the scale suppressing the quark and lepton operators should be the 
same.) 

One-loop WIMP-Quark Coupling Only: If no tree-level WIMP-quark opera- 
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Figure 4. Reach of the ILC dark matter searches, as a function of the dark matter mass M x 
and the thermally averaged annihilation cross section at freeze-out, (cr(xX ~* e + e~)v). The 
regions above the curves are accessible at the 3-sigma level. The value of (cr(xx ~* e + e~)t>) 
needed to obtain the correct relic density, in the absence of other annihilation channels, is 
also shown. 
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tors are present, the dominant contribution to WIMP-quark scattering appears at one 
loop, see Fig. 3. The LHC bounds are irrelevant in this situation due to the weakness 
of the interaction, so let us focus on direct detection. Since the momentum exchange 
in direct detection is well below Aqcd, the Z boson contribution can be neglected, and 
the photon couples to the nucleon as a whole, not to individual partons. Thus, only 
protons participate in this interaction. For simplicity, let us further assume that the 
WIMPs only couple to electrons (not /x or r), and that the coupling occurs through 
one of the operators listed in (2.2). With these assumptions, it is straightforward to 
evaluate the cross section of elastic WIMP-proton scattering, a p , for given M x and A, 
see Ref. [6]. Note that the interaction in Fig. 3 is unavoidable, and adding tree-level 
WIMP-quark couplings can only increase the cross section unless there are accidental 
cancellations. Thus, our a p can be thought of as an approximate model-independent 
lower bound on the direct detection cross section for given M x and A. 

In Fig. 5, we plot the contours of a v corresponding to the ILC reach for various 
run scenarios. Of the four operators considered in this paper, only the vector and 
t-channel scalar couplings give non-zero contribution to the diagram in Fig. 3. The 
other two operators would produce much smaller direct detection cross sections for the 
same A, and are not included in the plots. To compare with the published experimental 
bounds on WIMP-nucleon cross sections, one needs to rescale those bounds to take into 
account that in our model only protons, and not neutrons, couple to WIMPs, and the 
cross section depends on the momentum exchange. Including these effects, the current 
XENON100 bound rules out o p down to about 10 -43 cm 2 for masses above 50 GeV or 
so, with sensitivity decreasing rapidly for lower masses [6]. This leaves a large region 
of allowed parameter space which can be explored by the ILC. 

4 Conclusions and Outlook 

In this paper, we analyzed the sensitivity of the proposed high-energy e + e _ collider, 
the ILC, for a signal of direct WIMP pair production, accompanied by a photon. We 
employed the operator formalism to describe the WIMP-electron coupling, giving the 
analysis a high degree of model-independence. We found that the ILC-250, the "Higgs 
factory" stage currently under active discussion, can extend the reach in terms of A 
(the scale where the WIMP-electron interaction is induced) by a factor of 2.5 — 3 
compared to the current best bounds from LEP-2. Higher energy ILC upgrades can 
yield further spectacular improvements in sensitivity, especially if a beam polarization 
option is available. In all cases, the ILC reach allows to probe WIMP annihilation 
cross sections well below the value required for the thermal relic density to match 
the observed value. In other words, the ILC can discover WIMPs even if the process 
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Figure 5. Direct detection (WIMP-proton elastic scattering) cross section induced by the 
one-loop diagram in Fig. 3, for values of A at the outer edge of the ILC experimental sensitivity. 



XX e + e~ only contributes a tiny fraction of the total annihilation rate XX ~~ ^ SM 
in the early universe. We also discussed comparison with direct detection experiments 
and collider searches at the LHC While no firm model-independent conclusions could 
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be drawn, we argued that the ILC will be able to explore large regions of parameter 
space consistent with present bounds. 

If the 7 + tfi signal is discovered, the next task will be to measure the WIMP mass, 
the scale A, and to characterize the helicity structure of the operator coupling WIMPs 
to electrons. It would be interesting to study quantitatively the potential of the ILC 
to make these measurements. (See Refs. [3], [7] for similar studies in slightly different 
theoretical frameworks.) We leave this issue for future work. 
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Note added: 

While the writing of this paper was being completed, Ref. [12] appeared which 
overlaps with our work. 



A Analytic Formulae for Differential Cross Sections 

In this Appendix, we list analytic expressions for double-differential cross sections of the 
process e + e _ — » xXli m photon energy and its angle 9 with respect to the electron 
beam, for each of the four WIMP-electron interaction operators listed in Eq. (2.2). For 
convenience, we define 

- 2£l /.= ^. (A.1) 



\/S 

21 



Kinematically accessible range of z is [0, 1 — 4/x ]. The formulas below give the double- 
differential cross section, , < ? a , for polarized electrons and positrons: for example, 
"LR" corresponds to the e^e^, etc. For the case of vector interaction, we obtain: 



LR — RL : ^ ± -L_ ; l^i-^ (1 - , + V) [4(1 - ,) + * ( 1 + cos* *)] ; 
RR = LL : 0. (A.2) 
The case of axial-vector interaction yields 

RR — LL : (^^y 1 (' - ') Wl - *) + *V + »* ^ : 

LR = RL : 0. (A.3) 
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For the "s-channel" scalar interaction, we find: 



RR = LL : 



as 1 



(^)V.> ( * 



8tt 2 A 4 z sin 2 6 
LR = RL : ; 



+ - 2 ) - 



(A.4) 



while for the "t-channel" scalar interaction case, 



RR = LL : 



a 



1 



'1 - z - 4/i : 



2 r 



LR = RL : 



1927r 2 A 4 zsin 2 ^ \/ (1 - 2) 3 
-2/i 2 (cos 2 + 4(1 - z) - z(3z 2 -6z + 4)) 
as 1 



1 - z) ((2 - ^)(3z 2 - 6z + 4) - cos 2 9) 



96tt 2 A 4 z sin 2 6 



1-z-A^ 
l-z 



(l-z + 2/i 2 ) [4(1 -z) + z\l + cos 2 9)] . 

(A.5) 



We have checked that the Ward identities (vanishing of amplitudes upon replacement 
of the photon polarization vector by its momentum) work for all cases. 

We denote the polarization of the electron and positron beams, respectively, by P_ 
and P + . In our notation, — 1 < P± < 1; P_ = +1 corresponds to pure and P + = +1 
corresponds to pure e£, while P± = corresponds to unpolarized beams. With this 
notation, the cross section is given by 



l-P-\ fl-P+\ fl + P-\ fl + P+ 



+ 



2 

1-P_ 



2 

1 + P^ 



0"LL + 



2 

1 + P- 



2 

1-P, 



0"RR • 



(A.6) 
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